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THE  THEORY  AND  PRACTICE  CP  THE  h-p  VERSION 
OF  FINITE  ELEMENT  METHOD 


Ben  Qi  Guo*  Ivo  Babuska** 
Institute  of  Physical  Science  (  Technology 
and  Department  of  Mathematics 
University  of  Maryland 
College  Park,  MD  20742 
U.S.A. 


Abstract ;  The  p  and  h-p  version  is  a  new 
development  in  finite  element  method  in  recent  years. 

)  This  paper  is  addressing  some  theoretical  advances  and 
presents  numerical  illustrations,  t - — —  - r  ' 


For  a  survey  of  the  state  of  the  art  of  the 
and  h-p  versions  we  refer  also  to  (2). 


Finite  Element  Method  and  the  Approximation 


1.  Introduction 


Let  B(u,v)  be  a  bilinear  form  defined  on  H^*H2> 


There  are  three  versions  of  finite  element  method. 
The  classical  h-version  achieves  the  accuracy  by  re¬ 
fining  the  mesh  while  using  low  degrees  p  of  ele¬ 
ments,  p  »  1,2  usually.  The  p-version  keep  the  mesh 
fixed  and  the  accuracy  is  achieved  by  increasing  the 
degree  p.  The  h-p  version  properly  combines  both 
approaches . A 


where  are  reflexive  Banach  spaces  equipped  with 


norm  H  *|l  ^  and  H  "Uj*  resPectively  •  l*1  further 


F  t  H*,  i.e.,  F  be  the  linear  functional  on  Hj - 


We  seek  the  solution  uQ  c  such  that 


-  The  h-p  version  is  the  new  development  of  ^finite 
element  method.  It  was  first  addressed  by  Babuska  and 
Dorr  (4].  The  further  analysis  and  computation  for  two 
dimensional  problems  were  made  by  Guo,  Babuskap  In — ^3- 
vhere  the  exponentiaT  F5W~6f  convergence  was  proved. 

The  one  dimensional  analysis  was  given  by  Gui,  Babuska, 

— (181^-  The  improvement  of  the  results  for  curvilinear 
boundary  and  curvilinear  elements  was  made  by  Babuska, 
Cub.  Irr  fgJi  '  The  problem  with  non-hpmogeneous  Dirichlet 
data  was  studied  by  Babuska,  Guo.in  (10]  ^  The  h-p 

version  with  C^-elementa-f^r  the  problem  of  2m  order 
was  discussed  by  Guo.in  (191.  The  feedback  agg  adap¬ 
tive  approach  was  developed  by  Gui,  Babuska -HDK and 
Babuska,  Rank«^liy_.__ 


B(u0>v)  «  F(v)  for  all  v  e 


If  the  bilinear  form  B(u,v)  is  continuous  and 
satisfies  the  well-known  inf-sup  condition  (see  (31). 
then  the  problem  (2.1)  has  a  unique  solution. 


Let  now  c  h^,  S^  c  H^.  The  approximation 


problem  is  to  find  the  finite  element  solution  u  < 

S1 


such  that 


B(ug  ,v)  ■  F(v)  »  B(Uq,v)  for  all  v  c  S?.  (2.2) 


The  p  and  h-p  version  for  two  dimensional 
problems  were  implemented  in  the  commercial  code  PROBE 
by  Noetic  Tech.,  St.  Louis.  See  (27,28].  The  commer¬ 
cial  program  FIESTA  (Istituto  Sperimentale  Hodelli  e 
Struttre)  has  limited  p  and  h-p  capabilities  in 
three  dimensions.  The  p  and  h-p  versions  of  fi¬ 
nite  element  method  in  three  dimensions  is  being  devel¬ 
oped  by  Noetic  Tech,  and  by  Aeronautical  Research 
Institute  of  Sweden. 


For  sake  of  simplicity  we  assume  that  the  bilinear 
form  B(u,v)  satisfies  the  inf-sup  condition  on 
Sj^Sj.  Then  uQ  exists,  is  unique,  and 


l,us  "Vi  5  c(s1,s2)z(u0.h1.s1) 


Z(u0,H1>S1)  =  inf  Hu0-w||. 

wcS^ 


For  details,  see  (1,31.  We  shall  assume  that 


The  practical  effectivity  of  the  p  and  h-p 
version  is  closely  related  to  the  problems  in  structu¬ 
ral  mechanics  which  are  the  problems  of  elliptic  par¬ 
tial  equations  with  piecewise  analytic  data  (as  do¬ 
mains,  boundary  conditions).  The  analysis  of  the 
interior  regularity  was  given  by  Morrey  (231 ,  and  the 
behavior  of  the  solution  in  the  neighborhood  of  corners 
and  edges  of  the  domain  was  given  by  Kondrat'ev, 

Oleynik  (21,22]  and  Grisvard  (16,19].  The  characteri¬ 
zation  of  the  regularity  of  the  solution  on  nonsmooth 
domains  in  the  frame  of  countably  nonned  spaces  was 
ginven  in  the  series  of  papers  by  Babuska,  Guo  in  (5,6, 
7,8)  and  (9). 


C(SrS2)  5  D 


where  D  is  independent  of  S^  and  Sj,  hence  the 


norm  of  the  error  e  *  ug  -uQ  is  completely  governed 


by  ZiUg.Hj.Sj).  The  accuracy  of  finite  element  solu¬ 


tion  is  actually  an  approximation  problem.  0bvi6us)y 


the  solutions  u^  is  not  known,  but  it  will  be 


assumed  that  they  belong  to  a  compact  solution  set 
1C  <■  .  Hence  we  are  interested  in  Zfu^.H^.S^)  for 


any  u^  c  K.  The  precise  characterization  of  the  set 


The  advances  in  the  p-version  were  discussed  in 
(261 .  The  computational  aspects  were  addressed  in 
(25). 


K  which  is  made  as  small  as  possible  is  crucial  for 
the  most  accurate  approximation. 


3.  Characterization  of  the  Solution  for  Elliptic 
Problem  with  Piecewise  Analytic  Data 


*  Supported  by  the  National  Science  Foundation  under 
Grant  DMS-8S-16191 . 


••Supported  by  the  Office  of  Naval  Research  under 
Contract  N000014-8S-k-0169. 


As  indicated  in  the  previous  section,  selection 
of  the  h,p,h-p  versions  and  the  performance  of  three 
versions  strongly  depend  on  the  solution  set  K.  Let 
us  characterize  this  set  for  elliptic  problems  with 
piecewise  analytic  (input)  data. 


»0PY 

.  INSPECTED  , 


WHS 


'  Let  Q  e  ■  be  a  bounded  domain  (curvilinear 
polygon)  shown  in  Figure  3.1  with  vertices  ,  12 

i  S  M.  the  boundary  3ft  be  a  piecewise  analytic  curve 
H 

r  -  U  f,  where  r.’s  are  open  arc  with  endpoint  A. 
i-1 

and  Ai+J.  We  denote  the  internal  angle  by  0  < 

<  2«,  1  <  i  <  t 

A.  . 


g  e  Hk  V2(t,  (reap.  then  there  is 

G  £  hJ’So)  (reap.  b|(Q))  such  that  g|^  >  g.  We 
define 

l  %  i.  »  %  ..  *  inf  ffcB  .  ■  . 

k-1/2,1-1/2^,  c|  k.|Q) 


Theorem  3.1.  Let  ft  be  a  polygon,  f  t  Bg(0) ,  g*  t 
“Bi/2(r1),  g°  £  B2/2(r°),  0  <  $.  <  1.  8.  »  1  -  — 

p  6  1  1  u. 

(resp.  1  -  if  different  types  of  boundary  condi- 

i 

tion  are  imposed  on  I\  and  ri+1> •  1  2  i  S  M,  then 

the  problem  (3.1)  has  a  unique  solution  u  t  H2(ft)  and 
u  £  B2 (ft) .  For  proof,  see  [5). 


Fig.  3.1.  The  domain  with  piecewise  analytic  boundary. 

Let  r°  -  U  ?.  and  T1  -  r-r°  be  the  Dirichlet 
UD 

and  Neumann  boundary  respectively.  For  simplicity  we 
consider  only  the  problem 


If  ft  is  a  curvilinear  polygon  we  introduce 

C2(ft>  -  {u  e  Hg' 2 (ft) ,  |d°u(x) |  <  Cdk‘1(k-1) ! («k(x))_1, 

V  a,  |a|  ■  k  i  1,  C  and  d  independent  of  k ) .  Then 
we  have 


♦  u  *  f 

in 

(3.1a) 

Theorem  3 

0 

u  *  g 

on 

r° 

(3.1b) 

r.*s  are 

l 

3u  1 

aTT  ’  « 

on 

r1 

(3.1c) 

BB/2l<rl: 

«  4  (r1) 

0. 
.9  1 

L  t  H1  (f  ) , 

1  i 

i  *  D. 

1  "  2S71  ’ 

g°  is  continuous  on  r°,  the  problem  (3.1)  has  unique 
solution  (weak  sense)  ug  c  H1 (ft) .  No  matter  how 

smooth  f.g^g1  are,  the  singularity  appears  at  the 
corners  of  the  domain.  Hence  the  standard  Sobolev 
spaces  are  not  a  powerful  tool  for  this  type  of  prob¬ 
lem,  and  various  weighted  Sobolev  norms  were  intro¬ 
duced  . 

Let  6  =  (B^, . . ., B^)  be  an  M-tuple  of  real  num¬ 
bers  0  <  8^  <  1,  1  2  i  £  M.  For  any  integer  k  -  0 
vie  shall  write  B+k  =  (B.+k, . . .  ,8^+k) .  By  r^  (x)  we 

denote  the  Euclidean  distance  between  x  c  ft  and  the 
vertex  A.,  1  5  i  5  M.  We  denote  then  *..i.(x)  “ 

1  P+K 

M  3  +k 
n  r  1  (x). 

i-1 

Define  *or  k  >  *.  >  0,  Hk,fa)  -  (u  £  Ht_1(U), 

P 

#B+k-l0°u  *  (Q) ,  1  -  |a|  <  k)  (if  l  -  0,  the  con¬ 

dition  that  u  £  H*  * (ft)  is  absent)  and  bJ (ft)  - 

0 

(u  £  H^Sft),  <Cdk-l(k-l).,  k  -  t. 

t+1,  t+2 .  | ii |  -  k,  C  and  d  independent  of  k). 

As  usual,  we  denote  a  »  (a^.a^),  |a|  *  Oj+a.,,  ?  0, 

a  jl°lu 

i  »  1,2,  integers  and  D  u  « - -  u 

°1.  °2  °1  °2 
”X1  *X2  *1  X2 

Let  y  be  the  union  of  some  edges  of  ft.  The 

space  H*  *^2’*  1^2(fl)  (resp.  B*  ^/'2(y))  is  defined 
P  D 

as  the  trace  of  II* '*(0)  (resp.  bJ(P))  on  y.  If 

P  u 


1  -  t — )  ,  1  !  i  5  H.  Then  the  problem  (3.1)  has  a 
i 

unique  solution  u  £  Hg(0)  and  u  £  Cg(0).  For  proof, 
see  [9) . 

Remark  3.1.  Since  Bg(Q)  c  Cg(Q)  e  B2+t(fl)  for  arbi¬ 
trary  c  >  0  the  result  of  Theorem  3.2  is  weaker  than 
that  of  Theorem  3.1.  Nevertheless,  it  will  not  affect 
the  asymptotic  rate  of  convergence  for  the  h-p  ver¬ 
sion. 

Remark  3.2.  Theorems  3.1  and  3.2  are  also  valid  for 
generaly  strongly  elliptic  equation  and  system  with 
analytic  coefficients  satisfying  inf-sup  condition 
[5,6).  The  interface  problems  with  piecewise  analytic 
interfaces  and  the  eigenvalue  problems  have  the  same 
properties  too,  see  (8) .  The  solution  of  elliptic 
problem  of  2m  order  on  polygonal  domain  belongs  to 

b"+*(0)  |7) .  Hence  for  this  class  of  problems  includ- 

D 

ing  many  structural  problems  the  solution  set  K  « 

Bg (ft)  or  Cg(ft).  t  i.  2. 

The  definition  of  B«  2^2(y)  does  not  give  the 
0 

structure  of  the  space,  and  it  is  often  difficult  to 
verify  in  general  whether  g  belongs  to  this  space. 
Hence  further  characterisation  of  the  structure  of 

B*  '^2(y)  is  important  for  application. 

6 

.  .  .  ,  _ 1  _  ,  ,  .  . 


Let  I  -  (a.b)  c  ■  .  Analogously  as  before  we 

shall  define  the  spaces  H**'*(I)  and  Bl(I).  Let  o, 

y  y  : 

*  |x-a|,  0^  -  |x-b|  and  y  «  (Yj.Yj)  be  the  2-tuple 
of  real  numbers,  0  <  y  «  1,  i  -  i,2.  For  any  inte¬ 
ger  k  t  0  we  shall  write  y  «  k  1  (y.»k,y  «k),  anu 


(Yl*k,Y2«k) .  anu 


be  the  maximum  mk) 


H  Y^k 

dene!  c  T  •  II  p.  .  For  k  Z  t  2  0,  we  define 
T+k  i-l  1 

-  (u  «  Ht'l(I).  T  tu(,)  <  L2(I),  t  £  m  £k) 

(if  t  -  0,  the  condition  that  u  «  HX  X(I)  is  ab¬ 
sent)  and 

b‘(I)  -  (u  t  Hk,t(I).  llTY+)t.tu(k,|IL^(I)  £Cdk_t(k-t)! 

k  ■  1,1+1,...,  C  and  d  independent  of  k). 

Analogously  we  define  spaces  BX  (IV)  on  I\ 


with 

Yi 

instead 

of 

Y.  We  have  the  following  theo- 

rent. 

Theorem  .1 

ti- 

Let 

Q 

be  a  polygon. 

further 

(1) 

If 

0 

g 

is 

continuous  on 

„0  0  0, 

F  '  qi  "  9  If.  4 

■IjV'  1  4  D-  With  ri  ■  <Yi.l'Vi.2)'  0  4  Ti.t  4 

1/2,  t  -  1,2,  then  there  is  a  function  G°  c  B; (Q) 

such  that  G°|  Q  -  g°.^  and  6i  «  max  (tj  j)  * 

1/2  for  i,i-l  t  D,  Bi  »  1/2  ♦  2'  tor  1-1  4  D 

and  i  /  D,  0  <  B^  <  1  are  arbitrary  for  i,i-l  i  D. 


origin  and  0.  .,  h.  ,  and  h.  . 

i.l  i.j  -i.3 

minimum  of  length  of  edges  of  Q.  . ,  then  d.  . ,  h,  . , 

i,3  i,3  -i,) 


h.  .  satisfy 

A#  J 


n+2-j 


d.  .  <  on*x‘i. 

1.3 


*  ,d.  .  i  h.  <  h,  .  5  K  „d.  . 
1  i.j  i.j  2  1.1 


for  1  £  i  5  i(j),  1  <  j  £  n+1  and 
di.l  '  °* 

ic  on+1  <  h  <  h  <  (C  o° 

3  — i ,  1  i.j  ~  4 

for  1  £  i  £  i(l).  it  ,1  £  a  £  4  are  independent  of 
m 

i  and  j. 


(C2)  Let  M  -  (M.  .,  1  £  i  £  I(j),  1  £  j  £  n+1). 

1  »3 

H.  ^  is  a  one-to-one  napping  of  standard  square  S 

(resp.  standard  triangle  T)  onto  Q.  . .  Let  P  and 

1,3  i 

Y,  denote  the  vertex  and  side  of  Q.  . .  then  H  1 . (P  ) 

£  i»j  X , j  £ 

and  M.X.(y.)  are  the  vertex  and  side  of  S  (resp. 

l, 3  l 

T) ,  1  £  t  £  4  (resp.  1  £  t  £  3) .  Moreover,  if  M.  . 

1,3 

and  M  nap  5  (resp.  T)  onto  element  8.  .  and 

m, k  1,3 

Q  .  with  comon  side  y  •  A  A  ,  then 
W  , k  X  X  2 


(2)  If  gX  -  g1 | r  t  B°  (I*.),  with  y.  - 
i  Ti 

(Y-  .  ,Y-  _),  0  <  Y.  #  <  1/2,  i  €  D,  l  -  1,2,  then 

there  is  a  function  G1  <  BX(8)  such  that  GX|  ,  -  g1, 

6  rl 

and  Bi  -  max  (y1  *  1/2  for  i'i_1  *  D’  ^  " 

y.  .+1/2,  for  i  ID.  i-l  «  D,  0  <  8.  <  1  are  arbi- 
1  ,  X  1 

trary  for  i,i-l  c  D. 

Remark  3.3.  Theorem  3.3  also  holds  for  curvilinear 
polygon  with  piecewise  analytic  boundary  (see  [10]). 

Remark  3.4.  Theorem  3.3  allows  us  to  verify  the  boun¬ 
dary  conditions  mentioned  in  Theorems  3.1  and  3.2. 

4,  The  Mesh  and  Finite  Element  Spaces 

Mesh  design  is  very  crucial  to  the  accuracy  of 
method  and  depends  very  much  on  the  solution  set  K  » 

B2(S)  and  <^(Q).  We  assume  for  simplicity  that  0 

..a  a  polygon  contained  in  a  unit  disc  centered  at  ori¬ 
gin  which  coincides  with  the  vertex  A,  of  8,  and 
2  a  * 

K  -  B_(Q)  with  ■  r  ,  i.e+*,  assume  that  the  sin- 

O  D 

gularity  appears  only  at  one  vertex  of  Q. 


dist  (M?1.  (A),  mT1.(A.))  «  dist  (M_1.  (A),  m"1.  (A,)), 

l,j  !»}  *  BI,K  IM ,  R  £ 


for 


any  A  €  Y^,  £*1,2.  We  assume  each  side  y^ 


of 


Q.  .  is  analytic  curve,  12154  (resp.  IS  12  Z) , 
*■ ,  J 


hi.3*i.3,l<X> 

hi.3Vi.3.t(X) 


x  e  X  -  (0,1) 


and 


U(k)  I. 

1,  j 

where  C  and  L  are  independent  of  l, i.j.  Accord¬ 
ingly,  the  mapping  M.  .  of  S  (resp.  T)  onto  8.  . 

i,3  1.3 

is  analytic  on  S  (resp.  ¥)  and  can  be  extended  to 

S*  3  S.  Let  J.  .  be  the  Jacobian  of  M.  ..  We  shall 

1.3  i.l 

assume  that 


,.(k> 


I*'-;  |  £  ClSi 

1  ,  J  •  *■ 


C.h.  .  £  J.  .  £  C,h.  . 

1-1,3  1,3  2  i,J 

with  constants  C^.Cj  independent  of  i.j. 

Remark  4.1.  Figure  6.3  is  an  example  of  the  geometric 
mesh  Tor  the  problem  with  singularity  at  one  corner, 
but  the  mesh  can  be  analogously  generalized  for  prob¬ 
lems  with  singularity  at  every  corner. 


Mesh  typically  used  in  the  h-p  version  is  such 
that  domain  is  divided  into  several  layers  by  geometric 
progression,  the  jth  layer,  1  £  j  £  n+1  consists  of 
c'.  . ents  Q  ,  l£i£l(j).  In  addition  to  the  usual 

1.3 

conditions  in  the  theory  of  finite  element  method,  the 

mein  characterization  of  the  (geometric)  mesh  P1^  * 

(a.  .,  1  £  i  £  Kj),  1  £  j  £  n+1)  is  following) 

i,3 


(Cl) 

<  0  <  1, 


Let  mesh  factor  o  be  an  arbitrary  number, 

and  let  d  .  be  the  distance  between 
i.j 


Remark  4.2.  If  mesh  8n  contains  triangular  elements 
■  „  , -  o 

some  additional  assumptions  have  to  be  imposed.  In  the 
practice  these  assumptions  can  easily  be  satisfied,  see 
(9). 

Let  P  -  (p.  . ,  1  £  i  £  I ( j)  ,  1  £  j  -  n+1)  and 

i,3 

(?•  (q.  .,  1  £  i  £  Kj),  1  £  j  £  n+1)  be  the  degree 

1,3 

vector  with  integer  p.  .  and  q.  .20. 

i.j  1.3 

We  define  the  finite  element  spaces 


6.  numerical  Miatti 


SP'°((£l  -  (♦(♦<». yl  -  <*•»»»  *<»  <*.»»  « 

fl.  4((>M  is  the  polynomial  of  degree 

*  *  J 

5  j  in  C  *nd  of  dtgrM  '  j  An 
n) 

and 

sP.O.ifgOj  .  sp.9(an,  „  H^O), 
a  a 

Ip'0,1<fin)  •  Sp'3(fl")  n  H1(Q) 
o  o 

where  H^di)  »  (u  c  hNo),  u|  »  0). 

r 

By  N  we  denote  dim  (sP’®’1  (Q^) ) ,  the  number  of 
degree  of  freedom. 


We  will  present  feme  numerical  results  for  the 
solution  of  a  plane  strain  elasticity  problem.  We 
selected  the  model  of  crack  .panel  loaded  by  traction 
that  the  exact  solution  is  the  first  (symmetric)  and 
second  (antisymmetric)  node  of  stress  intensity  factor 
solution.  This  problem  was  selected  bee  use  it  charac¬ 
terises  the  usual  difficulties  of  engineering  computa¬ 
tion.  Due  to  the  syssaetry  and  antisymmetry  we  need 
only  to  solve  the  problem  in  the  upper  half  of  the 
panel  shown  in  Figure  6.1.  The  solution  has  singular 
behavior  at  the  tip  of  the  crack,  i.e. ,  the  displace¬ 
ment  U  -  (u,v)  has  the  expression  (r1^2^)*), 

rl^262(#))  near  the  origin.  Obviously  u,v  {  H2(0) 

and  u,v  c  Bg(0)  for  6  >  1/2. 


5.  Basic  Approximation  Theorems  of  the  h-p  Version 


The  energy  of  U  is  defined  as 


We  will  list  some  basic  approximation  results  in 
the  case  that  H2«  H1(Q),K  -  B2(Q)  or  C2(Q)  and 

S,  ■  S.  ■  sP,0'1<nn> ,  i.e,  we  seek  the  estimates  of 
12  o 

2(u,H1(0),  SP,0'1((l"))  for  u  <  1C. 
o 

Theorem  5.1.  Let  Q  be  a  polygon  and  u  t  B2(fl>  n 

(Q) ,  then  for  any  a  e  (0,1),  P  -  Q,  vj  <  p^  <  pn, 
0  i  v  <  u  <  «  and  p  ,  £  1.  We  have 

Z(u,H1((l),  SP'Q,1(an)>  <  Ce-bN  /  (5.1) 

0 

where  b  and  C  are  independent  of  N  * 

dim(SP'®'1 (o”) ) ,  the  number  of  degree  of  freedom.  For 

proof,  see  (9). 

Theorem  5.2.  If  u  e  C2(0)  n  0  1*  a  curvilinear 

polygon,  the  boundary  of  domain  is  piecewise  analytic, 
then  the  result  of  the  previous  theorem  holds.  For 
proof,  see  [9], 

Remark  5.1.  Mesh  factor  a  can  be  any  number  t  (0,1) 
the  computation  shows  that  a  ■  0.15  is  the  optimal 

value.  In  (18J  it  has  been  proved  that  o  *  (VJ-l)2 
•  0.17  is  the  optimal  mesh  factor  in  one  dimensional 
setting.  The  value  a  -  0.15  in  two  dimensional  prob¬ 
lems  reflects  the  fact  the  solutions  in  the  neighbor¬ 
hood  have  essentially  one  dimensional  character. 

Remark  S.2.  If  g?  «  g°|.  t  b}  (F-,),  i  i  t  with  y, 

1  1 i  'i  1 

»  (B^-l/2,  B^+j-l/2)  are  non -homogeneous  Dirichlet 
boundary  condition,  the  theorems  above  hold  provided 
g°  is  properly  projected  on  the  trace  of  finite  ele¬ 
ment  space  SP'0'1(0on). 


G(U) 


2 (l-2v) (1+v) 


/{(l-v)((|*)2  ♦  (|^)2) 


»y 


d/2-v,((g  +  ^)2)dxdy  -  ||0H^ 


where  E  and  v  are  the  roung’s  modulus  of  elasti¬ 
city  and  Poisson  ratio.  The  error  e  ■  0  -  is 

measured  in  energy  norm,  and  by  (5.1) 


HeHE  <  Ce 


-bN 


1/3 


The  relative  error  is  defined  as 

Be||E.R  “  l«lE/*u»E,<100%- 


t 


Figure  6.1.  The  crack  panel. 

The  computation  of  the  p  and  h-p  version  were 
made  by  program  PROBE  (27).  The  computation  of  the 
h-version  with  p  •  1  was  made  by  the  adaptive  pro¬ 
gram  FEARS  developed  at  University  of  Maryland  (24). 

We  will  compare  the  perfonsance  of  the  three  versions 
of  finite  element  method. 

Meshes  1  !  i  !  6  which  are  refined  near  the 

tip  by  geometric  progression  with  factor  o  •  0.15 
are  shown  in  Figure  6.2. 


Remark  5.3.  For  problems  of  order  2m,  the  theorems 
hold  when  geometric  mesh  contains  only  parallelogram 
and  triangular  elements.  For  details,  see  (22). 


Table  6.1.  Relationship  between  |«IE  .  M.n.p.b  and 
C  for  the  h-p  version.  The  symaietric  problem 


MESH  A, 


MESH  A. 


Figure  6.2.  Mesh 


The  table  6.1  shows  the  relationship  between 

is  the  number  of  layers. 


N.n.p  and  HeHE  R  where  n 


p  is  the  element  degree  and  M  is  the  number  of 
degree  of  freedom .  The  relationship  is  plotted  in 

Ln  ||e|  scale  and  shown  in  Figure  6.3.  The 

curve  characterising  the  convergence  of  the  h-p  ver¬ 
sion  is  the  envelope  of  the  six  curves  of  the  p- 


version  on  the  Mesh  A^,  1  £  i  5  6  and  is  nearly  a 


straight  line.  The  slope  of  the  line  is  b  in  (6.1) 
and  is  numerically  0.67. 


Figure  6.3.  Relative  error  in  energy  norm  vs.  number 
of  degree  of  freedom.  The  symmetric  problem  (E  »  1 . 
V  »  0.3)  on  Mesh  A  .  1  5  i  *  6,  o  *  0.15,  p  ■  n. 


In  Figure  6.4  we  show  the  dependence  of  the  error 
on  o.  We  see  the  best  value  of  o  is  close  to 

(^2-1) 2  ■  0.17  which  is  the  theoretical  optimal  value 
in  one  diawnsion. 


Figure  6.4.  Dependence  of  Relative  error  of  the  h-p 
version  in  energy  on  the  mesh  factor  o.  The  anti¬ 
symmetric  problem  (E-l,v-0.3)  on  Mesh  A^ ,  15156 


Figure  6.5  shows  that  the  h-p  version  in  insen¬ 
sitive  to  change  of  Poisson  ratio.  The  slope  of  the 
curves  of  the  h-p  version  for  v  *  0.3  and  0.49 
are  almost  the  same.  The  locking  problem  never 
occured. 

N - - 

25  CO  200  400  €00800000 


Figure  6.5.  Insensitivity  of  Relative  error  of  the 
h-p  version  in  enerqv  norm  to  change  of  Poisson  ratio 
The  antisymmetric  problem  (E=l.v=0.J)  on  Mesh  A 


Table  6.2.  Estiaatad  error  of  the  h-p  version.  The  symmetric  probl 


v  -  0.3)  on  Mesh 


1  «  1  <  6,  a  m  0.15.  p 


n 

Mr 

Mr 

Mr..*- 

Mr..*- 

(Mr-M<VMc*-  j 

1 

2.9596E-I 

2.9662E-I 

60.83 

60.92 

0.2189  3 

2 

9.8774E-2 

9.851  IE-2 

20.28 

20.23 

-0.2669  1 

j 

3.7033E-2 

3.705SE-2 

7.606 

7.611 

0.0606  ' 

0.0926  J 

4 

I.2489E-2 

I.2500E-2 

2.565 

2.567 

5 

4  3359E-3* 

4.3691  E3 

0.891* 

0  897 

06689  : 

—  * 

In  Figure  6.6  we  compare  the  performance  of  the 

h.p  and  h-p  versions  in  In  |e|_  xln  N  scale.  We 

E,R 

see  that  the  accuracy  0 . 5~1 . 0»  is  very  expensive  and 
probably  is  not  achievable  at  all  for  the  p-version 
and  h  version  with  p  »  1.  The  h-p  version  allows 
us  to  use  a  relatively  very  small  number  of  elements  to 
obtain  high  accuracy. 


25  30  PO  200  400  800 


h-VERSUN  I 

*csh - 1-2 


I  KT*1'! 


h-PVERSCN  ' 
_  MESH  A-Z 


3  4  3  6  7 

LN  N  - 

Figure  6.6.  Relative  error  in  energy  norm  vs.  number 
of  degree  of  freedom  for  the  h.p, h-p  versions.  The 
symmetric  problem  on  which  mesh  factor  o  -  0.15. 

Since  the  exponential  rate  of  convergence  can  be 
achieved  for  low  p  and  n  as  well  as  high  p  and 
n,  we  can  use  the  computational  results  from  three 
successive  computations  to  obtain  estimation  of  error 
in  energy  norm.  Table  6.2  shows  that  the  estimated 
error  is  very  reliable  with  relative  error  less  than 
1%.  For  the  formulation  of  estimated  error,  see  (20). 


7.  Conclusions 

From  theoretic  analysis  and  computation  we  can 
conclude  the  following: 

(1)  The  asymptotic  theory  reflects  accurately  the 
computational  practice  in  the  entire  range  of  engineer¬ 
ing  accuracy.  The  exponential  rate  of  convergence  is 
achieved  in  computation  in  industrial  practice. 


(2)  The  optimal  geometric  factor  of 
ment  is  close  to  0.15. 


esh  refine- 


(!)  The  performance  of  the  p  and  h-p  versions 
is  not  influenced  when  the  Poisson  ratio  »•  Si  0.5 
(i.e.,  when  the  material  is  nearly  incomprehensible).  | 

(4)  Industrial  experience  with  the  method  (by 
program  PROBE)  indicates  high  effectivity  and  advan¬ 
tages  of  the  h-p  version,  see  115). 


(S)  Preliminary  computation  and  theoretical  ana¬ 
lysis  show  that  in  the  three  dimensions  the  p  and 
h-p  have  superior  qualities  in  practical  computation 
of  problems  in  structural  mechanics. 
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